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We report on a measurement of the branching fractions for B~ — > D^*'''^ 'k' l~ vt and — > 
O'^*^'^ ■k'^ l~ vi with 275 million BB events collected at the T(4S) resonance with the Belle de- 
tector at KEKB. Events are tagged by fully reconstructing one of the B mesons in hadronic 
modes. We obtain B{B' D+n-rue) = (0.54 ± 0.07(stat) ± 0.07(syst) ± 0.06(BR)) x 10"^ 
B{B- D*+7v-£-ui) = (0.67±0.11(stat) ±0.09(syst) ±0.03(BR))xlO-^ B{BO ^ D^'-K+rve) = 
(0.33 ± 0.06(stat) ± 0.06(syst) ± 0.03(BR)) x 10"^ B{B^ D*'^7v+rve) = (0.65 ± 0.12(stat) ± 
0.08(syst) ±0.05(BR))xl0"^, where the third error comes from the error on B ^_D^*'>£-ue decays. 
Contributions from _B" — » D*^ t.~vi decays are excluded in the measurement of _B" — » 0*^11^ fvi. 



PACS numbers: 13.20.He 

Semileptonic decays play a prominent role in the 
study of B meson properties. The total semileptonic 
branching fraction has been precisely determined to be 
(10.73 ± 0.28) X 10-2 Q. While B Di'Di and D*i-i?i 
account for 70% of this total, other contributions are 
not yet well understood. The most promising candi- 
dates include resonant and non-resonant D^*^tt in the fi- 
nal state. Both ALEPH ^] and DELPHI have studied 
B — > D^*^'K£~Pi decays (where both 77=*= and 7r° modes are 
included). Assuming that semileptonic B decays other 
than B D£~vi and B D*£~D( are of the form 
B -> D^*\£-Dt, they find: 

= B{B^ DTTrv,) + B{B~^ D*tt£-v,) 
= (2.16 ± 0.30 ±0.30)xl0"2 (ALEPH), 
= (3.40 ± 0.52 ±0.32)xl0"2 (DELPHI). 

The former result suggests that there is a significant un- 
known contribution to the semileptonic branching frac- 
tion, while the latter shows no such deficit. More precise 
measurements are therefore desired to resolve this dis- 
crepancy and clarify the difference between the rate for 
the inclusive semileptonic decay and the sum of the rates 
for the exclusive modes. Improvement in the knowledge 
of the B — > D^*^'k£~i>i branching fractions will also help 
to reduce systematic uncertainties in measurements of 
Cabibbo-Kobayashi-Maskawa elements such as \Vcb\ and 
\Vu,\ 0. 

In this paper, we present measurements of the branch- 
ing fractions for B — > D''*^7r£^Pe decays. Inclusion of 
charge conjugate decays is implied throughout the pa- 
per. The analysis is based on data collected with the 
Belle detector |5j at the KEKB e+e" asymmetric col- 
lider We use a 253 ib'^ data sample at the T(4S') 
resonance (-^/s ~ 10.58 GeV), corresponding to a sam- 
ple of 275 million BB pairs. The selection of hadronic 
events is described elsewhere 0- An additional 28 fb"^ 
data sample taken at a center-of-mass energy 60 MeV be- 
low the T(4S') resonance is also used to study continuum 
e+e" ~* qq {q — u, d, s, c) events. 



The Belle detector is a large-solid-angle spectrome- 
ter with a 1.5 T magnetic field provided by a super- 
conducting solenoid coil. Charged particles are measured 
using hits in a silicon vertex detector (SVD) and a 50- 
layer central drift chamber (CDC). Photons are detected 
in an electromagnetic calorimeter (ECL) comprised of 
CsI(Tl) crystals. Kaon identification is performed by 
combining the responses from an array of aerogel thresh- 
old Cerenkov counters (ACC), a barrel- like arrangement 
of time-of-fiight scintillation counters (TOF), and dE/dx 
measurements in the CDC. A K/^k likelihood ratio Pk/ttj 
ranging from (likely to be a pion) to 1 (likely to be a 
kaon), is formed. With the requirement Pk/tt > 0-6, 
the kaon efficiency is approximately 88% and the aver- 
age pion mis-identification rate is about 8%. Electron 
identification is based on a combination of dE/dx in 
CDC, the response of ACC, shower shape in ECL, and 
the ratio of energy deposit in ECL to the momentum 
measured by the tracking system. Muon identification is 
performed using resistive counters interleaved in the iron 
yoke, located outside the coil. The lepton identification 
efficiencies are about 90% for both electrons and muons 
in the momentum region above 1.2 GeV/c, where leptons 
from the prompt B decays dominate. The hadron mis- 
identification rate is less than 0.5% for electrons and 2% 
for muons in the same momentum region. 

We use a GEANT-based Monte Carlo (MC) simula- 
tion to model the response of the detector and deter- 
mine its acceptance i3 ^ D^*^£~ve and D^*Jtt£~vi 
events are modeled using the EvtCen program 9]. We 
use an HQET-based mod el [loll for B D^*H~vi and 
the Goity-Robcrts model [lllffor B D^*'>Tr£-Vi. B 
D**{~^ D^*\)£^vt is also simulated using the ISGW 
model (l^ to evaluate the model dependence. 

To suppress the high combinatorial background ex- 
pected in the reconstruction of final states including a 
neutrino, we fully reconstruct one of the B mesons, re- 
ferred to hereafter as the tag. This allows us to separate 
particles created in the tag decay from those used in re- 
constructing the semileptonic decay, which we call signal. 
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It also provides a measurement of the momentum of the 
signal B meson, thus greatly improving the resolution on 
the missing momentum. 

The tag is fully reconstructed in the following modes: 
B+ ^ ^(*)07r+, D(*')°a+, D(*'>°Di*'>+ , and 

^ 7^W-a+ Z?W~Di*)+. 

candidates are reconstructed in seven modes: 
DO K+TT~, K+tt-tt", K+tt+tt+tt-, K^tt^, K^tt+tt-, 
ii'gTr+TT^Tr^, K^K^. D~ candidates are reconstructed 
in six modes: K'^tt^tt^, +7r+7r^7r°, K^tt^, 

if^7r~7r°, Kg7T^TT~7T~ , K^TT^ . candidates are 

reconstructed in D+ Kp<+, K+K-it+. The D can- 
didates are required to have an invariant mass to^i within 
±4 — 5cr of the nominal D mass value depending on 
the mode. D* mesons are reconstructed in the Z?*+ 
D"n+/D+TT", D*o ^ D^TTO/DO-f and D*+ D+-f de- 
cays. D* candidates from modes that include a pion are 
required to have a mass difference Am = ruDTr — i^d 
within ±5 MeV/c^ of its nominal value. For the de- 
cays with a photon, we require that the mass difference 
Am = mrt^ ~ ran be within ±20 MeV/c^ of the nominal 
value. 

The selection of the tag candidates is based on Afbc = 

V^bcam ~ Pb ^'^d AE = Eb - -Ebeam, whcrC ^boam = 

^/s/2 ~ 5.29 GeV, and pb and Eb are the momentum 
and energy of the reconstructed B in the T{AS) rest 
frame, respectively. The background from jet-like con- 
tinuum events is suppressed on the basis of event topol- 
ogy: we require the normalized second Fox- Wolfram mo- 
ment to be smaller than 0.5, and jcos^thl < 0.8, 
where ^th is the angle between the thrust axis of the 
B candidate and that of the remaining tracks in an 
event. The latter requirement is not applied to _B+ — > 
DOTT+,D*°{^D°Tr°)TT+ and B°-^D*-{-^ D°tt-)tt+ de- 
cays, where this background is smaller. 

The signal region for tag candidates is defined as 
A/bc > 5.27 GeV/c^ and -80 MeV < A£; < 60 MeV. If 
an event has multiple B candidates, we choose for each 
of B'^ and B^ the candidate with the smallest based 
on the deviations from the nominal values of AE', tod, 
and Am if applicable. 

Figure ^ shows the distribution in A/bc of and 
B^ candidates in the AE signal region. The number of 
tagged events is obtained by fitting the distribution with 
empirical functions: a Crystal Ball function for signal 
and an ARGUS function for background The fits 

yield (4.26 ± 0.17) x 10^ B+ , with a purity of 55%, and 
(2.72±0.11) X 10^ B°, with a purity of 50%. These yields 
include cross- feeds between _B+ and B^^ these peak in 
A/bc and are not well separated by fitting. From the MC 
simulation, we estimate the fraction of B^ [B^) events in 
the reconstruction of B+ {B") to be 0.095 (0.090). The 
corrected yields of tags are: iVtag = (3.88 ± 0.20) x 10^ 
for B+ and iVtag = (2.46 ± 0.12) x 10^ for B°. 

On the signal side, we reconstruct the following modes: 




5.24 5.25 5.26 5.27 5.28 5.29 5.3 5.24 5.25 5.26 5.27 5.28 5.29 5.3 

{GeVIc) {GeVIc) 

FIG. 1: A/bc distributions for fully reconstructed hadronic (a) 
and (b) B'^ decays. The solid curve is the sum of the fitted 
signal and background components, the dashed curve is the 
fitted background component. The signal region is indicated 
by solid arrows. 

B- D^*^"e-Dj, D(*)+7r-£~j>^, and B^ -> D'^*^+£~Di, 
B -> D^*H~vi decays are reconstructed 
as control samples to determine experimental resolutions, 
evaluate systematic errors, and normalize the results. All 
modes are reconstructed with all remaining particles after 
the full reconstruction of the tag. Flavor combinations of 
the two B mesons are restricted to B+B', B^B", B^B" 
and B^B^ . We do not require opposite flavor for neutral 
B mesons because of the effect of mixing. 

The selection of D^*^ mesons is performed in the same 
manner as for the tag. We require that lepton candidates 
be identified as electrons or muons and have momentum 
greater than 0.6 GeV/c in the laboratory frame. The 
electron and muon with largest momentum are selected 
as lepton candidates in each event. The pion candidates 
are selected with a loose criterion, Pk/tt < 0.9. We se- 
lect the B meson candidate formed with the best /?*^*) 
candidate (smallest of deviations from the nom- 
inal values of mo, and Am^j if applicable), and require 
that no charged track be left unused. Additionally, in the 
case of B° DOtt+I-Di and D*0{-^DOttO ,DO-f)TT+l~Vi 
reconstructions, we require A/(/I'°7r+) > 2.1 GeV /<? to 
veto events from B^ D*~^£~i>i. 

The semileptonic decay is identified by the missing 
mass squared, A/^j^^ = (£^bcam - - Ei- f - 

(Pstag + PD(')(7r) + Pi f, where £;boam is the beam en- 
ergy, /?£)(.) (jr) and El- are the D'-*\t:) and lepton en- 
ergies, and P^(,)(jr) and Pg- are the corresponding 3- 
momenta. Pst^ is the 3-momentum of the tag. All 
these variables are calculated in the T(4S') rest frame. 
For the signal, A/^^jg^ peaks around zero. The resolution 
on A/j^jgg ranges from 0.03 to 0.07 GeV^/c*, depending 
on the mode. In contrast, previous analyses (e.g. by 
the ARGUS collaboration [l6j), in which the Pstag term 
was neglected, led to a resolution of ^^0.5 GeV^/c'* on 
A/^jgg. This very good M^-^^ resolution allows reduction 
of the combinatorial background and to separate semilep- 
tonic decays that differ by only one pion or 7 in the final 
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FIG. 2: Af^i3^ distributions for (a) B~ -> D^Tvi, (b) 
B° D+Tvi, (c) B~ D'^rvi and (d) B° ~> D*+rue. 
Data are plotted as a histogram, and the fit result is overlaid 
as a dotted line. Signal (solid curve), other semileptonic de- 
cays (dash-dotted) and misidentified hadronic events (shaded 
histogram) are also shown. 

state. The distributions in Mj^^j^g for B D^*H~Vi and 
'vi are shown in Figs. 12 and |31 The signal events 
are clearly evident. 

We have investigated the backgrounds from BB and 
continuum events using MC simulation and off-resonance 
data. There are four main sources of background: 

1. Semileptonic B decays vifhere a pion or a photon is 
missed (e.g. B D*{Tr)£~De reconstructed as i? ^ 
D{Tr)£~Di if the soft tt*' or 7 from D* is missed). 
These distribute at high values of Af^j^g and can 
therefore be distinguished from signal. 

2. Semileptonic B decays where a random photon can- 
didate is used in reconstruction (i.e. B — ^ 
D°{tt)£~vi is reconstructed as i? ^ D*^{-k)1~ Di). 
These events have a lower value of M^j^^. 

3. Hadronic decays B D'^*^nh where one of the 
hadrons is misidentified as a lepton. These events 
peak near Afj^jj^g — and therefore need special 
care as described below. 

4. Random combinations. As the tag is fully recon- 
structed, this background is expected to be small. 
The distribution in M^j^^ is studied using events 
in the /\E side-band (0.1 GeV < A£; < 0.3 GeV) 
and off-resonance data. It is found to be consistent 
with a flat distribution in Af^j^g. 

The signal yields are obtained by a binned maxi- 
mum likelihood fit to the M'^^- distributions in the in- 



FIG. 3: A/^i,, distributions for (a) B" D+TTTue, (b) 
B" ^ D°TT+£-ue, (c) B- D*+-K-£-Vi and (d) B° 
D'^'tt'^ £" ui. Distributions are shown as described in Fig.|2| 

terval [-1.0,1.0] GeV^/c'^. The total fit function is: 
Fix ^ Ml,J = (7V,ig + 7Vbkg3)5(x) + 7VbkgiSi(x) -I- 
Nh]<ig2B2{x) + N\j-kg4B4{x), where iVgig is the number of 
signal events and the A^bkg(n) are the number of back- 
ground events corresponding to categories 1-4 listed 
above. The signal shape S{x) is a sum of two Gaussians 
and a single-sided exponential convolved with a Gaus- 
sian; the smeared exponential is necessary to describe 
the upper tail, which is due to energy loss, mainly by 
radiation from the electron. Bi{x) are normalized back- 
ground shapes. For Bi{x), we use a threshold function 
(oc exp(— /3a;)), plus an exponential function to de- 
scribe the lower tail component at Af^j^^ ~ GeV^/c"*. 
B2{x) is a smeared exponential function. B^i^x) is con- 
stant in the fitted interval. 

A'sig, A'bkgi and A^bkg4 are floated in the fit, while A^bkg2 
and A^bkg3 are determined from data and fixed. (Accord- 
ing to our background categorization, A^bkgi (A^bkg2) is 
non-zero only for B — > D[-K)£~h't {B — > D*'^{-k)£~ D^) de- 
cays.) The shape parameters of Bi{x) and B2{x) are 
fixed using MC simulations. The mean and width of 
S{x) are floated in the fit to i? ^ D^*'>£''p£, and the ob- 
served discrepancies with MC simulations are used to fix 
the mean and width of S{x) for B D'-*^n£-De. Other 
shape parameters of S{x) are determined from MC sim- 
ulations. 

To determine iVbkg2 we multiply the measured num- 
ber of i? ^ D^{n)£~Vi events by the probability 
of such events faking B — > D*^{-k)£~vi, as deter- 
mined in MC simulations. To determine the contri- 
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bution from misidentified hadronic decays (A^bkgs), we 
use a sample of events which are reconstructed as sig- 
nal but where the lepton candidates do not pass the 
identification requirement. We weight these events 
by lepton mis-identification rates binned in laboratory 
momentum, as determined from a control sample of 
D*+^D°{-^K-tt+)tt+ events. (In the case of 5° 
D(*)+7r~£^t?£, we also consider the possibility that the 
prompt pion is misidentified as a lepton.) iVbkgs is then 
determined by a fit to the M^^j^g distributions of this spe- 
cial sample. 

The results for B -> D^*h-vi and B -> D^*^7Te~Di 
are shown in Fig. ^ and listed in Table HI |Il| The 
statistical significance of the signal yield is defined as 
E = ■\/2 ln(— Lo/£max), where L^ax and Lq denote the 
maximum likelihood value and likelihood value obtained 
assuming zero signal events, respectively. Signals are ob- 
served in all B — > D'^*^'k1^vi modes with a statistical 
significance of more than 7. 



Mode 


N ■ 


^bkg3 




s(io-") 


PDG S(10-") 


DTve 


1049 ± 44 


9±4 


5.70 


2.37 ±0.10 


2.15 ±0.22 




1419 ±59 


13 ±4 


3.02 


6.06 ±0.25 


6.5 ±0.5 




467 ± 26 


4±2 


4.15 


2.14 ±0.12 


2.14 ±0.22 


D*+rvi 


476 ± 25 


3±2 


2.05 


4.70 ± 0.24 


5.44 ±0.23 



TABLE I: Signal yields, number of misidentified hadronic 
events, efficiencies, and raw branching fractions for B 
D^*^£~ui, compared with their world average values. Sys- 
tematic errors are not included. 



Mode 




-^bkg3 


6(%) 


E 




142.1 ± 16.5 


5.7 ±3.2 


3.41 


12.4 




62.5 ±9.7 


2.6 ± 1.8 


1.38 


9.9 




72.0 ± 12.4 


1.1 ± 1.7 


4.06 


8.2 




62.7 ± 11.6 


1.7 ± 1.5 


2.09 


7.2 



TABLE IL Signal yields, number of misidentified hadronic 
events, efficiencies, and statistical significance for B 

The raw branching fractions are calculated as S — 
iVsig/(2eA^tag), where A'sig is the measured number of 
B —> D'^*\-K)i~ vi events and A'tag is the number of se- 
lected tags. The efficiency e is defined as e = e^^isnai ^ 

signal / generic i generie • r n -l -l j x 

^frec /^frec ' whcrc e|j.^^ IS fuUy reconstructed tag- 
ging efficiency and Cf^f" ^' is its efficiency in the case of one 

B decays into the signal mode. The factor Cfrf" '^V^froT'^"^ 
accounts for the difference of B tagging efficiency be- 
tween the signal decay modes and generic B decays and 
is estimated to be around 1.05 depending on the mode. 
The obtained raw branching fractions in normalization 
modes, B D^*H- 'vi are shown in Table These agree 
with the world average values quoted in Ref . P] . 



To calculate the branching fractions for B 

'1^1, we first find the ratio R of the raw branch- 
ing fractions for each of these decays to that of the 
B D^*'>£~Di decay mode with the same D'-*^ charge, 
then multiply R by the world average values of B{B ^ 
D''*H-vi) [J. The results are shown in Table Iml 

Systematic errors in the measurement of the branch- 
ing fractions for B D^*^Tr£~Pe are associated with the 
uncertainties in the signal yields, tag yields, and recon- 
struction efficiencies. Most of the systematic errors re- 
lated to the reconstruction of the D^*-' and the lepton, as 
well as the branching fraction of the D^*^ decays, cancel 
out in the ratio R. 

Systematic errors in the signal yields are estimated 
from uncertainties in the signal shape, background shape 
and number of misidentified hadronic events. For the 
signal shape, the mean and width in B — > D^*^'Ki~ i>i 
are shifted by their respective errors obtained from the 
control sample, and determined to be 2—4%, depending 
on the mode. For the background shape, we estimated 
the uncertainty by using a different shape. For example, 
the main background components in _B — > D'Kt~vi are 
changed from B 0*^1' Ui to B ^ D**{-^ D*TT)rDg. 
The uncertainty is 12% for B'^ D'^-rr^t^vi due to larger 
background events, while it is smaller than 4% for other 
decay modes. For the number of misidentified hadronic 
events, we estimated the uncertainty to be 2—3% by vary- 
ing A^bkgs by its error. 

The systematic error due to the uncertainty in the 
number of tags and amount of flavor cross-feed have to 
be considered as well, since the branching fraction calcu- 
lation involves ratios of different B meson flavors. The 
former is estimated to be 4% by varying the background 
shapes used in the Mbc fit, and the latter is estimated to 
be 3% by varying the fraction of cross- feed by its error. 

The effect of uncertainty on reconstruction efficiencies 
mostly cancels since we take the ratio to a sample that 
differs from the signal by only one pion. We, therefore, 
assign a total of 1% error due to tracking efficiency (based 
on a study of partially reconstructed D* decays) and 
particle identification (based on a study of kinematically 
selected D*+ D^{-^ K^n'^)n+ decays). The uncer- 
tainty due to finite MC statistics used to model the signal 
is 2-3%. 

To estimate the uncertainty in efficiency due to the 
modeling in the MC simulation, we compared B — > 
D**{^ D(*\)l-vt and B D^*K£-i?i. The difference 
of 10% was assigned as the error due to this effect. 

The total uncertainty is the quadratic sum of all above 
contributions, and amounts to 13—20%, depending on 
the mode. In the measurement of the absolute branch- 
ing fractions, we quote an additional 4—10% systematic 
error due to the error on branching fractions of the nor- 
malization modes B D''*^I^vi. 

As a cross-check an independent analysis was per- 
formed with tighter requirements on Z?(*' reconstruc- 
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tion, Ai5, and Mbc of the tag and a slightly different 
fitting procedure and efficiency calculation. Results of 
both analyses are consistent. 



Mode 


R 


S(10"") 


D'+Ti-e.-vi 


0.25 ±0.03 ±0.03 
0.12 ±0.02 ±0.02 
0.15 ±0.03 ±0.03 
0.10 ±0.02 ±0.01 


0.54 ±0.07 ±0.07 ±0.06 
0.67 ±0.11 ±0.09 ±0.03 
0.33 ± 0.06 ± 0.06 ± 0.03 
0.65 ±0.12 ±0.08 ±0.05 



TABLE III: Ratios and branching fractions for B —* 
D^*''TT£~i'e decays. The first error is statistical, the second 
is systematic. For B, the third error is due to the branching 
fraction uncertainties of B — > _D'*^^~i/f. 
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We compute the total branching fractions of i? — > 
'vg assuming isospin symmetry, B{B — s- 
D(*\^frVt) = ^B{B D^*'>Tr^£-De), to estimate the 
branching fractions of D^*'>7r'^ final states. We obtain 

l3DMne-p,{B~) = (1.81 ± 0.20 ±0.20)xl0-^ 
BoM^i-^^iB") = (1.47 ± 0.20 ±0.17)xlO-2, 

where the first error is statistical and the second is sys- 
tematic. Our measurements are consistent with the 
ALEPH result and significantly smaller than that of 
DELPHI. This clearly shows, as suggested by ALEPH's 
result, that the missing branching fraction in semilep- 
tonic B decays is not fully covered by these excited states. 
Further searches must be carried out in i? ^ D^*^-K'Kll~i>i 
modes. 

A study of the mass structure of D'-*^tt in B ^ 
D^*^Tri~i?i decays is needed to understand the decay me- 
chanics of i? — > D**£~i>i, which would provide crucial 
tests of HQET and QCD sum rules 0. Belle has re- 
cently observed aU four D** states {D[,Dn, Di and D2) 
in the hadronic B decay B — > D**tt [l^]- The corre- 
sponding semileptonic decay B D**i~Pe, however, 
has been observed only for B'^ —^ £'i(— > D*^TT^)£~i>i 
by CLEO . These contributions could be clarified by 
the method with fully reconstructed tags with a larger 
data set. 

In conclusion, we have measured the branching frac- 
tions of B- D(*)+7r-€-z?f, and B° D^*^°tt+ l-vi. 
These decay modes have been clearly observed in a clean 
environment thanks to full reconstruction tagging, and 
the direct measurement of these branching fractions was 
achieved for the first time. 
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